Quantum Dynamics with Mean Field Interactions: 

a New Approach 

Laszlo Erdos* , Benjamin Schlein^ 

66 '■ 

O Institute of Mathematics, University of Munich, 

^ ■ Theresienstr. 39, D-80333 Munich, Germany 

April 23, 2008 

<\ 

m ■ 

(N ; 

Dedicated to Jiirg Frohlich on the occasion of his 60th birthday, 
^ \ with admiration and gratitude 

^ (— | Abstract 

> . 

We propose a new approach for the study of the time evolution of a factorized iV-particle 
bosonic wave function with respect to a mean-field dynamics with a bounded interaction potential. 
'— The new technique, which is based on the control of the growth of the correlations among the 
| | particles, leads to quantitative bounds on the difference between the many-particle Schrodinger 

dynamics and the one-particle nonlinear Hartree dynamics. In particular the one-particle density 
■^j- ' matrix associated with the solution to the iV-particle Schrodinger equation is shown to converge to 

r***> . the projection onto the one-dimensional subspace spanned by the solution to the Hartree equation 

with a speed of convergence of order 1/N for all fixed times. 

O ; 1 Introduction 

oo . 

We consider a system of N interacting bosons in v dimensions described on the Hilbert space TLn = 
\ L^(M. Nu , dx\, . . . dxjv)) the subspace of L 2 (W Nu , dx\ . . . dx^) consisting of permutation symmetric 
^ ■ wave functions. Here the variables xi, . . . ,xn £ ^ u refer to the positions of the N particles. The 
Hamiltonian is given by 



N N 

H N = J2 ("Ax, + U( Xj )) + - V( Xl - X j) (1.1) 

j=l i<j 

and, under suitable conditions on the potentials U and V, acts as a self-adjoint operator on H^. 
The coupling constant 1/N in front of the interaction characterizes mean field models; it guarantees 
that the kinetic and the potential part of the Hamiltonian are typically of the same order. 

The dynamics of the system is governed by the iV-particle Schrodinger equation 

id t ip N , t = H N ip Njt Tp N ,t = e- iH » t ip N . (1.2) 
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In particular we will be interested in the time evolution of factorized initial data iPn{x\, . . . , xjv) = 
1X7=1 f( x j)- ^ turns out that, because of the mean field nature of the interaction, for large N, 
factorization is approximatively preserved by the evolution (|1.2j) . In a sense to be made precise, we 
have ipN,t — l ff N ) f° r a suitable <pt G L 2 (R"). If factorization were indeed preserved, it is simple to 
check that (ft must be determined by the solution to the self-consistent nonlinear Hartree equation 

id t ipt = -Aipt + (V * \<pt\ Z )<Pt, (1-3) 

with initial data ipt=o = ( P- 

A more precise analysis shows that the solution ipN,t can be approximated, in the limit N — > 00, 
by products of solutions to the nonlinear Hartree equation fjl .3j) only if we are interested in quantities 
depending non-trivially only on a finite number of particles. In other words, we cannot expect the 
L 2 -difference between ipN,t and <pf N to converge to zero. What can be expected and, for a large class 
of potentials, has also been rigorously proven, is that for an arbitrary fixed /c-particle observable 
(for example, a compact operator on L 2 (R fcy )), and for every fixed time t £ M, 

& N)t , (J« ® 1^-*))^,*) - (iff N , ( J (fc) l^ k ^f N ) = (vf\ ) (1.4) 

as N — > 00, where (ft is the solution to (jl.3p with initial data (ft=o = f - The convergence f)l .4[) can 
be interpreted as the convergence of the marginal densities. Recall that, for a given TV-particle wave- 
function TpN, we define the density matrix 7^ = \^n)(iPn\ as the orthogonal projection onto ipjy, 

(k) 

and, for k = 1,...,N, we define the A:-particle marginal density 7^ associated with ip^ by taking 

(k) 

the partial trace of 77V over the last N — k particles. In other words, 7)^ is defined for k = 1, . . . , N, 
as a non-negative trace class operator on L 2 (M, uk ) with kernel given by 

7^ ) (x fc ;x' fc ) = J dx N - k j N (x k ,x N _ k ;x.' k ,x N ^ k ) = J dx N - k if) N (x k ,x N - k )ip N (x k ,x N _ k ). (1.5) 

Here and henceforth we use the notation X& = (xi, . . . , x k ), x' k = (x[, . . . x' k ), x^v-fc = (x k +i, ■ ■ ■ , %n)- 
We will also use the shorthand notation x = x^r = (xi, . . . , xjv)- Note that we choose the convention 
TV 7^ = 1 for all N > 1, and 1 < k < N. 

(k) 

For k = 1, ... ,N denote by the A:-particle marginal density associated with the solution 
to the iV-particle Schrodinger equation (jl.2p with factorized initial data ipx = cp® N . In terms of 
marginal densities, (jl.4p is equivalent to the convergence 

y { £l^\<p t )(<Pt\ 9k asiV^oo (1.6) 

for all fixed t € R and k > 1, with respect to the weak* topology defined on the space of trace 
class operators. It is simple to check that the weak* convergence in this case is equivalent to the 
convergence in the trace norm topology, because the limit is a rank one projection. 

The first rigorous results of the form (jl.4p or (jl.6p were obtained by Hepp in [15j . and then 
extended by Ginibre and Velo in |14j . These works are based on the analysis of the time evolution 
of coherent states (in a second quantized representation of the system). The use of coherent states 
(characterized by a non-fixed number of particles) makes it possible to isolate the main part of the 
evolution (described by the Hartree equation) and to study the fluctuations around it. With this 
method, it was recently proved in [20] that (jl.6p can be improved (for example for the case k = 1) 
to the quantitative estimate 

TN,t- \<Pt){<ft\ < const (1.7) 
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for factorized initial data (while for coherent initial data, the results of [20J establish a better bound, 
proportional to 1/iV) and for interaction potentials V such that V 2 (x) < C(l — A x ). 

A different approach was proposed by Spohn in [21J. The main idea of the technique introduced in 

• (k) 

|21| consists in characterizing the limit of the marginal densities 7^ t through the evolution equations 

that they satisfy. For finite N, the dynamics of the marginal densities {TjvlltLi ^ s governed by the 
so-called BBGKY hierarchy 



k ^ 



\t 

7=1 i<j 

k ( L8 ) 

, N — k \ r ( fc+1 )- 

+ N 2^ Tr * i + 1 r ^ j ~ Xk+l ^ ry kt 

3=1 

where TVfc +1 denotes the partial trace over the (fc + l)-th particle. Using the a-priori bound Tr = 

1, an abstract argument shows that the sequence {^NUk^l is compact with respect to a suitable 
weak topology. The first main step towards a proof of (|1.6|) in [21J is then to show that an arbitrary 
limit point {7«^j' ) j}fc>i of the sequence of marginals {7jvi}fcLi is a solution to the infinite hierarchy 

k k 

idtltlt = E [" A ^S] + E Tr *+i [ V ( x i - ^+i)>^J 1} ] (1-9) 

3=1 j=± 

which is obtained from (jl.8p by formally letting N — » 00. The second main step to show (|1.6p consists 
in proving the uniqueness of the solution to (jl.9p . Since it is simple to verify that the ansazt 7^ = 
\(ft){tpt\® k is indeed a solution to (|1.9p if tp t satisfies the Hartree equation (jl.3p . the compactness of 
the sequence {^\}^ =1 leads to (11.61) . This method, which was first used in [21] to prove (]1.6p for 
bounded interaction, was recently extended to singular potential (see [HI EIEIIIIIIIEIIIIEIEI E]). 
Note that this argument does not provide an effective estimate for the speed of convergence in (jl.6p . 



Recently, a new proof of (jl.6p has been presented by Frohlich et al. in [12]. This method provides 
convergence estimates uniformly in Planck's constant up to an exponentially small remainder, and 
therefore it can be used to combine the mean field- with the semiclassical limit (note that, in the 
semiclassical limit the Hartree equation converges to a Vlasov equation). In [13] it has been shown 
that the limit (jl.6p can be interpreted as a Egorov-type theorem, in the sense that the time evolution 
of observables commutes with their Wick quantization, up to corrections which vanish in the mean- 
field limit. 

In the present work we propose yet a different approach. With respect to the techniques mentioned 
above, our approach has the advantage of being very simple and of giving a quantitative bounds of 
the order 1/N for the speed of convergence in (jl.6p . The estimates we obtain here are therefore 
better, in their ^-dependence, than the bounds fll-7j) derived, for factorized initial data, in |20j . 
However, our analysis is restricted to a smaller class of interaction potential V; in fact, we need to 
assume that the operator norm [|V|| is finite (the operator norm of the multiplication operator V 
equals the L°°-norm of the function V(x)), and that it has a Fourier transform V € L 1 (M I/ ). The 
main result of this paper is the following theorem. 

Theorem 1.1. Fix an arbitrary dimension v > 1. Suppose that V € L°°(M. U ) is such that \\V\\i < 00. 
Denote byyffi t the k-particle density associated with the solution tpN t t = e~ tHNt ipN of the Schrddinger 
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equation with initial data ipN = 9?® N , for some cp G ff 1 (M 3 ), with \\<p\\2 = 1- Then we have 



Tr 



^ ^ v --/>■ I'll f e 8|[V||t _ ]_ 



/or every t £ M, N > 1 and 1 < k < N . Here we defined \y = 1 + (|| V||i/||V||) and ^ is the solution 
to the nonlinear Hartree equation 

idttpt = -A<ft + (V * \(pt\ 2 )pt , 

with ipt=o = <f- 

In order to prove Theorem 1 1.1 1 we control of the growth of the correlations among the particles (at 
time t = 0, the wave function is factorized, and therefore there are no correlations). Our proof has 
been inspired by the techniques used in [19] to prove Lieb- Robinson inequalities for anharmonic lattice 
systems (see also the proof of Lieb- Robinson bounds for quantum spin systems given in [TT 1 [T8 | [T6] ) . 
To prove Lieb-Robinson bounds, one has to control the growth of the correlations between spatially 
separated observables. To prove Theorem 1 1.1^ on the other hand, we have to control the growth of 
correlations between observables acting on different particles. 

In Section we prove that if A and B are two operators acting on different particles, then the 
commutator [A,e Ht Be~ iNt ] remains of order 1/N for all fixed times (the commutator is zero at 
time t = 0); see Theorem 12. 11 We apply this result to prove that, if denote the marginal density 
associated with the evolution ipN,t = e~ lHNt i^N, for a factorized initial data ipjf = tp® N , then the 
difference 7j^ + — Jpp) ®7jvl rema i ns of order 1/N for all fixed times (it is zero at time t = 0), when 
tested against product observables; see Proposition ^. 21 Note that both Theorem 1 2 . 1 1 and Proposition 
12.21 only require bounded potential to hold (the condition ||V||i < 00 is not required here). Finally, 
in Section [3l we apply the result of Proposition 12.21 to show Theorem 11.11 

Notation. For k > 1, denote by Bk the algebra of bounded operators over L 2 (M. ku ). For an 
operator A € B^, and for integers 1 < i\ < i% < ■ ■ ■ < i^ < N, we denote by A^ 1 '---' 1 ^ the operator 
on L 2 (M Nv ) acting as A over the particles i±, . . . ,ik and as the identity over all other particles. We 
will use the notation \\A\\ to indicate the operator norm of A (where A is meant as an operator from 
an I? space to itself). In particular ||V|| denotes the operator norm of the potential (interpreted as 
a multiplication operator) which of course equals the L°°-norm of the function V(x). 

2 Growth of Correlations 

Theorem 2.1. Suppose that V G L°°(W). Then, for every bounded m-particle observable A G B m , 
and every bounded n-particle observable B G B n , we have 

A [h,...,i m ]^ e iH N t B [j 1 ,...,j n ] e -iH N t < mn\\A\\ \\B\\ ^ e 4||y|||t| _ ^ ^ 2 _ 10 ^ 

for all integers 1 < i\ < ■ ■ ■ < i m < N, 1 < ji < ■ ■ ■ < j n < N, with {h,... , i m } n {ji, . . . ,j n } = 0- 
Proof. Without loss of generality, we assume that t > 0. For n, m G N, we define the quantity 

II r^[ti,...,i m ] e iH N tg[j 1 ,...,j n ] e -iH N t~\ II 
f m ,n{t)= SUp ' ii. (2.11) 
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for arbitrary integers 1 < i\ < ■ ■ ■ < i m < N, and 1 < j\ < ■ ■ ■ < j n < N such that {i±, . . . , i m } f] 
{iii ■ ■ ■ >Jn} = 0- Note that, because of the permutation symmetry of Hn, the quantity f(t) is 
independent of the choice of the indices ie,je and 



fm,n(t) = SUp 1 

We define next the modified Hamiltonian 



^[ii+l,..,n+tn^ e iH N tg[l,...,n] e -iH N t~^ | 



=1 j>n 



With respect to the dynamics generated by H K N , particles 1, . . . ,n are decoupled from the rest of 



An) 



the system. In particular this implies that e jv '■■■ ,n 'e l n * is still an operator acting only on 
the degrees of freedom of particles 1, . . . , n, with norm equal to the norm of B. Therefore, we have 

j^[n+l,...,n+m]^ e iH N t e -iH^- ] t £>[!,. ..,n] e i.ffjy* ] t e ~iH N t 



fm,n{ t ) = SU P 

■n 



B\ 



(2.12) 



For given A € B m and B £ B n we define the time-dependent bounded operator acting on L 

g AB (t) = J \ln+^,---,n+m]^ e iH N t e -iH^h^[l,...,n] e iH l ^ ) t e -iH N t 



and we compute its time-derivative 



j9A,B(t) = \Al n+1 >~> n+m \e iH » t [i(H N - H ^) }e -iH^t B [i,...,n]^h 



-iHpft 



^[n+i,...,n+m] ^ i e iH Nt(£[ N _ H^)e~ iHNt , e iHNt e~ i ^ B^ 1, '''' n ^e i ^ e~ iHlvt 
ie lH ^(H N - HP)e- iHNt ,g A:B (t) 



+ 



e iH N t e -iH { ™ ) tg[l,...,n\ e iH]t l l >t e -iH N t ^ j^n+l,...,n+m\ ^ & iH N t _ H^)e 



(n)x -iH N t 



where, in the last step we used the Jacobi identity. Next,we define 

H (n) {t) = e lHN \H N - H ( ^)e- iHNt . 

It is simple to see that (t) generates a two-parameter group of unitary transformations (i, s) 
satisfying 

id t U {n \t, s) = n (n) {t)U {n \t, s), with U {n \s, s) = 1 for all s € R . 
Therefore, we obtain 
d 



dt 



= 14^(0, t) e iH Nt e -iH < £ ) t B [l,...,n] e iH^'t e -iH N t 



An). 



A [n+l,...,n+m}^ e iH N t( HN _ H ^ e ~iH N t 



(2.13) 



W (n) (i,0) 
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Integrating this identity from time to time t, using that gA,B (0) = and the definition of we 
find that 



9A, 



N 



e=i j>n 



JH N s g-i-H"^ 1 ' s £j[l,...,n] e «-H"^ l) s e -iH N s 



(2.14) 



A [n+l,... I n+m.] j ^-^(^ _ Xj )e~ iHNt 



U {n \s,t) 



Taking norms, we get 



\9A,B 



<—J2J2 ds W B W A^ n+l — n+m \e iHNS V{x l -x j )e 

£=l j>n^° 

Amnt\\A\\\\B\\\\V\ 



< 



N 



+ 2n\\B\\ / As 



o 



A [n+l,...,n+m]^ e iH N sy( Xi 



)e~ iH 



N s 



where the first term on the last line corresponds to the terms with j = n + 1, . . . ,n + m, while in the 
second term we used the permutation symmetry. From (|2,lip and (|2.12|) . it follows that 



f minify — 



Amnt \\V\ 
N 



+ 2n\\V\\ [ dsf m . 2 (s) 
Jo 



(2.15) 



Iterating this equation for k times, we obtain that 



4mn||V|| 4ranl|V| 



fc-1 



t + 



N 



+ 2n\\V\\U\W 



N 



\fc-i 



r=l 



dsi . . . 



d s t s v- 



(2.16) 



Sfc-l 



'0 J 

With the a-priori bound f m %{&) < 2, it follows that 



dSfc fm,2{s r 



mn \ -« 

Jm,n(t) — / , 



r=0 



(r + l)! + fc! S iV !< 



1 +n 



fc! 



for all fc 6 N. Since the l.h.s. is independent of fc, we obtain (|2.10p . 
The following proposition is a useful consequence of Theorem 12.11 



□ 



Proposition 2.2. Suppose that V G L°°(R U ). Denote ipN,t = e _iffjv ViV solution to the N- 
particle Schrddinger equation with factorized initial data ipN = <£® N , for some <p € L 2 (M 3 ) with 
\\<p\\2 = 1- Then, for any A £ B m , -B € £> n; we /iaue 



(^v.t,^ 1 '---'^^' 1 '--^ 

mnllAII 111?! 



< 



e 8||V||* _ j 



(2.17) 



/or arbitrary integers 1 < n < • • • < i m < N, 1 < j 1 < ■ ■ ■ < j n < N with {i x , . . . , i m }n{ji,. . . ,j n } 
0. In particular, if 7^ denotes the k-particle marginal associated with ipN,t, we have 



Tr A 



[l,...,m] a R[m+ l,..., m+ n]^ ^m+n) _ Jm) ^ 



< 



mn\\A\\ \\B\ 



N 



e *\W\\t_^ 
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Proof. Because of the permutation symmetry we have 

= (^, t ,^ [1 '-' m] i? [m+1 '-' m+n] ^t) - (^ lt! A[ 1 -- m V JV)t )(te, J B [m+1 '-' m+ri Viv ) t> • (2.18) 
We observe that 

^ N)t A [1 - M B^ m+1 - m+n H N +) 

= (lfj N , e iHNt A [l,..M e -iH N t e iH N t B [ m +l,..., m +n] e -iH N t^ N ^ 



{iP N ,e iHNt A^ 



.,m] -iH N t 



\ip)(<ff N e iHNt B^ m+1 ' 



.,m+n] —iH N t 



A' 



+ J2(^N,e iHNt A^-^e- iHNt [l^M^" 1 ) ® (1 - |^>M) ® 1®^') 
i=i 

X e iH N t J^[ m + 1 ,-, m + n ] e - iH Nt^ jN ^ _ 

Since -i/w = V 9 ®^ we S e t 

(V'jv,t,A[ 1 '-' m ]s[ m+1 '-' m+ "Vjv,i) 

= (^ ii ,A[ 1 '-' m Viv,t)(^ t ,5 [m+1 -- m+n Viv,t) 



X (1 - | ¥) )( ¥J |)b1 )e iHjvt B [m+l,...,m+n] e -iJf N t 



m+n 



j=m+l 



N 



j=m+n+l 



+ Yj [e i%4 i [1 "" H e- i *' ) (1 - |y>)(¥>|)0] ( |^>M® (i_1) ® 

X (1 - \tp){^\) [j] , e iH N t B [m+l,...,m+n] e -iH N t ^ \ ^ 



(2.19) 



(2.20) 



where we used the notation (1 — |y)(¥>|)^ = 1^ 1 - > <8> (1 — \(p)(<f\)® 1^ J ' ^ for the operator acting as 
the projection (1 — |y)(<^|) over the j-th particle, and as the identity over the other (TV — 1) particles. 
Note that, by definition ((1 - \ip)(ip\)^) 2 = 1 - \<p){<p\. From Theorem[2J] we obtain 



< 



ran||A||||£|| , 4 ||v|| t _ s mnP||||fl|| f 4 ||V||t _ n 2 



N 



N 



mn\\A\\\\B\\ mu _ 
N 1 ' 



which completes the proof of the proposition. 



□ 
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3 Derivation of the Hartree equation 

Using the bounds on the correlations obtained in the previous section, we can now state and prove 
our main result. 

Proof of Theorem From the BBGKY hierarchy, we obtain the integral equation 



7g=" W (*)l^r-^E/ d S U( k \t-s) V( Xl - Xj ),^ 



N 



(N-k) 
N 



(3.21) 



3=1 



£ / dsW«(i - s) Tr k+1 V( Xj - x k+1 ),^ 



(fc+i) 



where we defined the free evolution 

i/(fc)( t ) 7 (fc) = e -»Ej =1 (-Ai+^))* 7 (fc) e iE*=i(-Aj+t/(* j ))* j 

and where T^k+i denotes the partial trace over the + l)-th particle. On the other hand, it is simple 
to check that, if (ft is a solution to the nonlinear Hartree equation 

id t (p t = -Apt + (V * \pt\ 2 ) ft 

then the orthogonal projection \(pt){(ft\® k satisfies the equation 

k r t 

10 



\<p t )(<p t \® k =UW(t)\ip)(<p\ m I dsUW(t-s)Tr k+1 Vixi-Xk+dMfal^V 

3=1 



Therefore, for an arbitrary ^-particle observable over L 2 (M. ku ), we obtain 



N 



i<3~° 



£/ AsTt (uW(s-t)jW) Vfa-x^-yW 



(k+l) 



+ i N^J ^TV [V( xj -x k+1 ).- x 
k rt 

(k+l) Jk) _ (1) 



3=1- Q 



£ / d, Tr (U^(s - t)jM) V( Xj - x k+1 ), 7^ - 7 « 7 £ 



£ / d S Tr (U^(s - t)jW) V(xj - x k+1 ), ( 7 £ - \<p s )(<Ps\® k ) ® \<P.)(<p. 



3=1 "° 
k ft 

iJ2 I dsTv (u^(s-t)J^) \v(x j -x k+1 ) n ik l®( 1 ( il-\ l p s ){i Ps 

3=1 







We bound the fourth and fifth term using that 

Tr (U^(s - t)jW) \v(xj - x k+1 ), (iP - \<p s )(if s \ m ) ® \<Ps)(<Pa\ 



<2||F||||J«||Tr tS-I^)^ 



(3.22) 



(3.23) 
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and that, analogously, 

Tr (U^(s - t)j( fc )) \v(x 3 - x k+1 )^\ ® ( 7 Sl " 

v j i \ /j (3 24) 

<2||y||||j( fc )||Tr 7 W -|^}(^| . 
As for the third term on the r.h.s. of (|3,22|) . we expand the potential in a Fourier integral. We obtain 



Tr (U^(s-t 



V( Xj - x k+1 ), h%t 1) - 7S ® 7^ 



.X 1 ) 



< 



< 



dg|V"(g)| Tr 



(1) 



< 



/" dg|F(g)| |Tr (u ik \s-t 

+ J dq\V(q)\ |Tr <"'"'' ' (^ (fc) (s-t 

2fc||^lii|[ t / (fc) li /j||v||. _ . 
iV ' 



(3.25) 



where we applied Proposition 12.21 Thus, from (|3,22p . we find that 
TrjW ( 7 $-|^r) 

<2A;||y||||j( fc )|| f dsTr 7 « - \<Ps}(<Ps\® k + 2k\\V\\\\J^\\ f dsTr ^\-\<p a )(<p a 
Jo ' jo ' 



(3.26) 



+ 



fc 2 IMIi||J (fc) ll / biivin A + 3^ 2 ll^lll|J (fc) l 



qv\\N 



N 



For a trace class operator A over L 2 (M !/fc ), let ||^4||i = Tr \A\ be the trace norm of A. If we denote 
by JC k the algebra of compact operators on L 2 (M. uk ), we have the variational characterization of the 
trace norm 

\\A\\! = Tr \A\ = sup |Tr J A\ . 
JeK k , || J||=i 

Therefore we find 

7$ - i^r* , < 2fcii^n Td, ( 7 <*> - i^x^r , + 



TS S - \<Pa)(<Ps 



+ 



k 2 Xy 

2N 



(3.27) 



,8||V||i 



1 



with Ay = 1 + 



i/\\V\\). For k = 1, we obtain 



7£J-|¥*><¥*I , < 4\\V\\ / d, 7 Si-b.><V. 



Ay 
! + 2iV " 



8||V||t 



Iterating this inequality n times, it follows that 



l ~ 2iV 



:,8||V'||f 



e -„. ,r _ i 

n-1 



\ " * /"i /-si rsm-i , 



,8||V||* 



(3.28) 





rsn-l 




■ / ds n 


/o 


JO 



72, Sn - l¥0<V««l 
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Since 



ft rs m -i 
jo jo 



ds m e 8 H y ll^ < - l 



and ll7S )Sn - \Ps n }{<P Sn \\\i < 2, we find 

Inserting this bound into the r.h.s. of (|3.27p . we conclude that 



(3.29) 



Ttyt ~ \<Pt){<Pt 
Iteration leads to 

t£ - \<Pt)(<Pt\ 



< 2k\\V\\ / ds 
i Jo 



Tjfl - \<Pa){<Ps 



k\v / 8||^||i _ -i 
l N 



k 2 Xy 

i - N 



< 



Mv\\t 



+ ^ L E( 2k wnr I d-i / d S2 



m=l 



+ (2fc||F||) n / dsi... 



/" 


ds 2 ■ ■ 


Jo 




ds n 


(k) 



ds m [e 



W\\s„ 



1 (3.30) 



which implies that 



7^-|^r| 1 <^e a * l|V|1 (e 8 '^-l) 



□ 
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